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ABSTRACT 


The  three  fundamental  boundary  value  problems  of  the  classi¬ 
cal  theory  of  linear  elasticity  have  been  formulated  In  a  manner 
that  permits  maximum  exploitation  of  numerical  techniques  and 
modern  high  speed  digital  computers.  Both  the  static  and  dynamic 
cases  have  been  Included  as  well  as  arbitrary  body  forces.  The 
extension  to  problems  Involving  more  than  one  media  has  been 
given  and  In  addition  no  restrictions  are  placed  on  the  boundary 


Introduction 

In  the  classical  theory  of  linear  elasticity  there  are  two 
classes  of  problems:  dynamic  and  static.  In  each  of  these  classes 
there  are  three  fundamental  boundary  value  problems  which  are 
labeled  the  first,  second  and  third  or  mixed  boundary  value 
problems  respectively.  The  first  problem  consists  In  determining 
the  displacements  In  an  elastic  body  when  the  stresses  are  speci¬ 
fied  on  the  boundary  and  the  second  problem  Is  the  determination 
of  the  displacements  throughout  an  elastic  body  when  the  displace¬ 
ments  are  given  on  the  boundary.  The  third  or  mixed  boundary 
value  problem  requires  the  determination  of  the  displacements 
In  an  elastic  body  If  the  stresses  are  prescribed  over  part  of  the 

boundary  and  the  displacements  are  given  over  the  remainder. 
(2),(3),(4)* 

Many  Important  problems  In  elasticity  may  only  be  solved 
by  employing  numerical  methods.  In  order  to  permit  maximum  use 
of  numerical  techniques  It  Is  desirable  that  all  three  boundary 
value  problems,  for  both  the  static  and  dynamic  cases,  be 
formulated  In  a  similar  manner.  Such  a  formulation  can  be 
achieved  by  employing  methods  analogous  to  those  used  In  potential 
theory.  In  this  paper  the  complete  formulation  of  these  problems 
will  be  given.  This  formulation  has  been  used  to  obtain  numerical 
results  to  some  specific  problems.  These  results  have  been 
given  In  reference  (5)  and  will  not  be  reproduced  here. 

♦Superscript  numbers  In  parenthesis  refer  to  references  In 
bibliography. 


Some  General  Considerations 


The  dynamic  equations  of  classical  linear  elasticity  may 
be  written  In  terms  of  the  displacement  vector, "u,  as 

+  fT  =  fUt« 

where  X  and  ^  are  the  Lame  constants,  Is  the  body  force 
vector  and  lower  case  subscript  letters  Indicate  partial  differ¬ 
entiation.  (6) 

u  Is  the  displacement  vector  with  components 
u,  V,  w  In  the  x,  y,  i  directions  respectively.  The  static 
equilibrium  equations  are  obtained  by  setting  the  Inertia 
force  terms  equal  to  zero  In  the  above  equations.  They  are 

+  ft  r  0  (2) 

Using  the  fact  that  a  vector  function  may  be  decomposed  Into  a 
solenoldal  and  an  Irrotatlonal  part,  the  displacement  vector 
may  be  written  as 

Iji  =  74  +  V  (3) 

where  =  0.  The  solenoldal  part,  V  >  corresponds 

to  a  rotation  and  the  Irrotatlonal  part,  V(J>  ,  to  a  dilatation. 

The  body  force,  F,  may  also  be  written  as 

where 
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and 


providing  that  P  Is  uniformly  continuous  and  vanishes  at  Infinity 
In  a  prescribed  manner. 

Substitution  of  Equation  (3)  Into  Equation  (l)  and  employing 
the  vector  Identities 


V  * 

< 

II 

o 

V- 

Vx^ 

=  o 

and 

V 

-  v('9 

r.Tp) 

-V‘y 

shows 

that  If 

and  satisfy  the 

Inhomogeneous 

wave  equations 

v'<i'  -I- 

^  °  -g  ‘f'tt 

and 

+ 

X  = 

where 

f  -  ^  - 

“i  ■  >1  +  2^ 

> 

X  = 

"A 

X  <-&« 

Ot  - 

> 

f 

then  the  displacement  vector  as  given  by  Equation  (3)  satisfies 
Equation  (l).  Similarly,  the  representation  (3)*  when^substltuted  lnt< 
Equation  (2)  shows  that  If  4>  and  'f'  satisfy  the  Poisson  equations 

V’tt*  +  f,  =0  and  Vf  +  X  =0  (5) 

then  the  displacement  vector,  Ti,  satisfies  the  static  equilibrium 
equations . 

Recalling  the  Integral  representation  of  the  solution  of 
the  Inhomogeneous  wave  equation,  the  functions  ^  and  may 
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be  written  as 


where  denotes  partial  differentiation  In  the  direction  of 

the  exterior  normal  to  (T  ,  r  is  the  distance  from  the  observa¬ 
tion  point  P  to  a  point  Pq  on  CT  and  (T  Is  the  closed  surface 
enclosing  the  volume  T  [^means  that  the  function  P(x,y,z,t) 

Is  to  be  evaluated  at  the  retarded  time  t'  =  t  -  r/c.  Equations 
(6)  and  (7)  are  Klrchoff's  solutions  of  the  Inhomogeneous  wave 
equations  In  terms  of  retarded  potentials.  The  corresponding 
Integral  representations  of  the  solutions  of  the  Poisson  equations, 
(5)  are 

and 
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If  body  forces  are  not  present  then  equations  (4)  and  (5)  sire 
replaced  by  the  homogeneous  wave  equations  and  Laplaces*  equa¬ 
tions  respectively  and  the  volume  Integrals  In  the  above  Integrals 
vanish. 

When  a  harmonic  time  dependence  of  the  form  Is 

assumed,  the  displacement  potentials  must  satisfy  the  following 
set  of  Inhomogeneous  Helmholtz  equations 

=  0  and  V'v  +  "Si  +  7  =0  (10) 


where 


^ 


and 


U'  -  CjD 

n»  “  “pf 

'-a 


The  Integral  representations  for  and  vy  are^^^ 

)  '  -S: 


and 


(11) 


(12) 
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As  before,  the  volume  Integrals  In  each  of  the  above  vanish  If 
there  are  no  body  forces  present. 

In  differential  geometry  It  la  shown  that  an  arbitrary 
surface  In  a  three-dimensional  space  may  be  characterized  by 
the  parametric  representation 


(13) 


where  and  (3  are  two  Independent  parameters.  If  the  lines  of 
curvature  at  a  point  of  the  surface  are  chosen  as  the  parametric 
representation  of  the  surface  then  a  local  orthogonal  coordinate 
system  may  be  constructed  at  the  point.  This  is  always  possible 
for  a  sufficiently  smooth  surface.  (10) 

The  coordinate  system 

consists  of  the  two  lines  of  curvature  and  the  normal  to  the 
tangent  plane  at  the  point.  The  unit  normal  Is  given  by 


y\  = 


(13) 


It  Is  assumed  that  the  surface  Is  sufficiently  smooth  so  that 
the  parametric  representation  as  given  by  Equation  (13)  possesses 
as  many  derivatives  as  may  be  required  In  the  subsequent  develop¬ 
ment  . 


In  the  next  section  It  will  be  shown  that  by  employing 
the  Integral  representations  of  the  functions  and  and 
Imposing  the  proper  boundary  conditions  that  a  set  of  Integral 
equations  may  be  derived  whose  solutions  are  the  desired  surface 
potentials  from  which  the  displacement  and  stress  fields  may  be 
calculated  by  Integration. 
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Formulation  of  the  Integral  Equations 

The  first  problem  to  be  considered  Is  the  second  fundamental 
boundary  value  problem  of  elasticity.  Let  D  denote  the  prescribed 
surface  displacement  vector  and  u  the  displacement 
vector  obtained  from  the  displacement  potentials .  On  CT  the 
normal  and  tangential  components  of  the  displacement  must  be 
equal  and  hence  on  (T  the  following  boundary  conditions  hold: 

=0-“^  and  =  D*n  .  (i4) 

In  terms  of  the  displacement  potentials  the  equations  are 

VcJ)  •  +  V  ^  (^5) 


and 


V<i> 


(16) 


In  a  previous  paper,  (5) 

the  following  vector  relationships  were 

derived: 


and 


Where  the  fact  that  Is  divergenceless  has  been  used  and 


H  =  sTes  -  F* 

with  E  - 

G  - 

F  =  X-T^  =  o 
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since  the  local  coordinate  system  Is  orthogonal.  Application 
of  these  vector  relationships  to  Equations  (I5)  and  (I6)  gives 


H4>«  +  =  HD-^  (17) 

and 

H  +  t.4>f  =  H'D  «7i  .  (18) 

Equations  (17)  and  (I8)  may  be  solved  for  the  surface  values  of 
and  In  terms  of  the  tangential  derivatives  of  and 
V  .  If  the  prescribed  surface  displacement  vector  Is  time 
dependent  the  values  of  and  are  substituted  Into 
Equations  (6)  and  (?)  to  give 


1  ^ 

an  \y\ 


-4-TT//niD-^nA)d-j? 


(19) 


and 
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y(Pjt) 


dt 


+  X.  *%  -%*  V.  ]  -  1}  dp 


-4ff'dt'D’‘'^l(T]d-dp 


(20) 

Allowing  the  observation  point  P  to  approach  the  boiindary  and 
Interpreting  the  Integrals  In  the  sense  of  their  limiting  values 
results  in  a  set  of  four  simultaneous  Integral  equations  for 
the  determination  of  the  surface  values  of  the  displacement 
potentials.  Solution  of  these  equations  for  the  surface  poten¬ 
tials  enables  the  field  values  to  be  calculated  by  Integration. 

This  Is  the  Integral  equation  formulation  of  the  second  funda¬ 
mental  boundary  value  problem  of  dynamic  elasticity.  If  the 
prescribed  surface  displacement  vector  Is  not  time  dependent  the 
problem  is  a  static  one  and  equation  (2)  Is  the  governing 
differential  equation.  Consequently,  the  surface  values  of 
and  %  obtained  from  Equations  (17)  and  (18)  are  substituted  Into 
Equations  (8)  and  (9)  to  give  the  required  Integral  representations 
for  and  Vj/  .  By  Interpreting  these  Integrals  in  the  sense 
of  their  limiting  values  the  Integral  equation  formulation  of  the 
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second  fundamental  boundary  value  problem  of  static  elasticity 
Is  obtained.  If  a  harmonic  time  dependence  Is  assumed  for  both 
the  prescribed  displacement  vector  and  the  body  force,  the 
values  of  and  'Vx  are  substituted  Into  the  Integral 
representations  of  the  Inhomogeneous  Helmholtz  equations.  The 
limiting  values  of  the  resulting  Integrals  gives  a  set  of  simul¬ 
taneous  Integral  equations  which  are  equivalent  to  those  obtained 
from  the  Integral  equation  formulation  of  the  scattering  of  steady 
harmonic  elastic  waves  from  rigid  surfaces  of  arbitrary  shape. 

The  next  problem  to  be  considered  is  the  first  fundamental 
problem  of  elasticity.  Let  T  denote  the  stress  vector  prescribed 
on  the  boundauy  O'  and  the  stress  In  the  body  acting  on  a  plane 
with  direction  cosines  1,  m,  n.  Then 

S  :r  Xt  +  Yj  +  I'll  (21) 

where 


The  stress-strain  relations  In  usual  notatlon^^^)  are 

(r«  zXe  +  2/^,  r'*  =/*(u,  +  ir,) 

(r«*=yw(Ui  +  W,) 

(r"=y<(wv+u,) 


(T"  =X0 

a"  =Xe+2/iur, 
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(23) 


where  the  dilatation,  0  ,  Is  given  by 


0  =v-n 


(24) 


The  boundary  conditions  require  that  the  normal  and  tangential 
components  of  the  stress  at  the  boundary  must  be  equal,  l.e. 


and 


=  T.ft 


(25) 


on  O'  .  (26) 

Substitution  of  Equations  (21)  and  (24)  Into  Equation  (25)  gives 

•tv  aI  .  /F<ui _ 1  I  /r** 


+2a’^J(w  +2o^*TT^‘n  +  =T-^ 


(27) 


Relating  the  stresses  to  the  displacement  potentials  by  using 
Equations  (3)  and  (25)  and  substituting  Into  the  above  equation 
yields 

+  (c5-2cO(4>«G  +ct)„E) 

In  a  similar  manner  Equation  (26)  may  be  written  In  terms  of  the 
displacement  potentials  to  give 

-  [Tij.  -  t,  {%•%)  +  %  ■%) 

--Tf 


(29) 
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The  algebra  necessary  to  obtain  equations  (28)  and  (29)  la 
quite  lengthy  and  Is  given  In  a  previous  paper. (5)  if  the 
prescribed  surface  stress  Is  time  dependent.  Equation  (l)  Is  the 
governing  differential  equation  and  the  displacement  potentials 
satisfy  Equations  (4) .  Since  the  local  coordinate  system  Is 


orthogonal. 


Ve 


Urp.(vl) 

(30) 

and  Equations  (4)  may  be  combined  with  Equations  (28)  and  (29)  to 
eliminate  and  ^Y«>\  giving: 

[‘♦‘-la  (4) 

#  Kfp  ©  ^1]  "*■  ^ 

(4-.S  +4>^6)  +  ?  ^ 


(31) 


and 
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(32) 


The  above  four  equations  when  combined  with  the  limiting  values 
of  Equations  (6)  and  (7)  constitute  a  system  of  eight  simultaneous 
equations  for  the  determination  of  the  eight  surface  potentials 
.  cti,  and  .  If  the  stresses  specified  on  the 
boundary  are  time  Independent  the  problem  Is  a  static  problem. 
Consequently,  the  time  derivative  terms  are  not  present  In 
Equations  (31)  and  (32)  when  Equations  (5)  are  used  to  eliminate 

and  from  the  boundary  conditions.  These  four  equations, 

together  with  the  limiting  forms  of  Equations  (8)  and  (9),  consti¬ 
tute  a  system  of  eight  equations  for  the  required  surface 
potentials  and  their  normal  derivatives.  Similarly,  If  a 
steady  harmonic  time  dependence  Is  assumed,  the  terms  ^it  and 
M'tt  in  Equations  (31)  and  (32)  are  replaced  by  and  — 

respectively  when  the  second  normal  derivatives  are  eliminated 
from  the  boundary  conditions  by  substitution  from  Equations  (lO). 
Again,  a  set  of  determinate  equations  for  the  surface  potentials 
and  their  normal  derivatives  may  be  obtained  by  adjoining  the 
limiting  values  of  the  Integral  representations.  Equations  (ll)  and 
(12).  These  sets  of  eight  equations  consltute  the  desired 
Integral  equation  formulation  of  the  first  fundamental  boundary 
value  problem  of  elasticity. 

The  third  fundamental  problem  of  classical  elasticity  Is 
formulated  In  an  euialagous  manner.  The  procedure  Is  a  direct 
extension  of  the  previous  formulations.  For  the  dynamic 
problem  the  limiting  values  of  Equations  (6)  and  (7)  must  hold 
everywhere  on  the  boundary  while  for  the  static  problem  the 
limiting  values  of  Equations  (8)  and  (9)  must  hold.  The  boundary 
conditions  require  that  on  those  sections  of  the  boundary  where 
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the  surface  displacements  are  prescribed.  Equations  (l?)  and 
(l8)  are  to  be  satisfied  and  on  the  remaining  sections  of  (T 
where  the  surface  traction  Is  given.  Equations  (28)  and  (29) 
must  hold.  These  latter  equations  can  be  combined  as  before 
with  Equations  (4)  to  eliminate  and  for  the  dynamic 

problem  and  combined  with  Equations  (5)  to  eliminate  the  same 
variables  for  the  static  problem.  Consequently,  at  every  point 
of  the  boundary,  there  exists  eight  equations  for  the  determination 
of  the  eight  surface  potentials  ct>  and  .  If  the 

dynamic  problem  Is  reduced  to  a  harmonic  time  dependence,  the 
Helmholtz  Integral  representations  are  to  be  used.  The  author 
would  like  to  point  out  that  parts  of  those  sections  of  Reference 
(5)  referring  to  the  second  and  third  fundamental  problems  are  In 
error  and  should  be  modified  In  accordance  with  the  analysis 
presented  here. 

For  those  problems  In  which  the  body  force  Is  absent  the 
displacement  potentials  satisfy  homogeneous  equations.  Thus,  In 
the  above  formulations  of  the  three  fundamental  problems  the 
volume  Integrals  are  not  present  In  the  Integral  equations. 
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Further  Modifications  and  Considerations 


When  the  elastic  medium  Is  a  composite  of  two  or  more 
different  elastic  media  the  previous  formulations  must  be  modi¬ 
fied.  The  procedure  for  setting  up  the  required  equations  for 
two  adjacent  media  will  be  outlined.  The  extension  to  more 
than  two  media  may  then  be  readily  accomplished.  Denote  the 
media  by  medium  I  and  medium  II  and  let  terms  In  medium  I  be 
denoted  by  the  subscript  I  and  those  terms  In  medium  II  be  denoted 
by  the  subscript  II.  In  addition,  let  denote  that  part  of 

the  Interface  between  media  I  and  II.  Continuity  of  normal  and 
tangential  displacements  across  O^n  Implies 

(33) 

emd 

(U  +  %  (<t), 

Continuity  of  normal  and  tamgentlal  stresses  across  Implies 

+  ((t'.L  (C^-£Ci)s  .  a),,  (Ci-2C4)E 
+  ('kl),,  (^li'  ^  Ci)  G  +  (.<X"  ^  ^ 
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and 


^  H  [1 X c/„  { 

+2  H  [Tu(<l>„),^  -  tp  (4„)..]  -  (^p  [(tLi; 

-7<[(Uf-^p]  ^^[{%)„-%] 


(36) 


where 


r  *  -.  ^1 ^A-  Xn +^/^rr 

">  — ft  Cb^,  -  V.: 


fa 


Ch-  -j- 


f-\i  _  x'^4*' 

'“Su  ■  « 


'w 


The  terms  (4^^^ ,  ,  (^)wn  and  (%)iw  may  be  eliminated 

from  Equations  (35)  and  (36)  by  employing  the  displacement  potential 
equations  appropriate  for  the  problem  and  the  medium.  The 
resulting  eight  equations  when  combined  with  the  limiting  values 
of  the  Integral  representations  yields  sixteen  equations  for  the 
determination  of  the  surface  values  of  the  sixteen  functions 
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In  addition,  on  those  segments  of  the  boundary  of  either  medium 
where  the  displacements  or  stresses  are  specified,  the  equations 
previously  developed  must  hold.  If  the  parameters  ol  and  ^ 
are  so  chosen  as  to  be  the  arc  lengths  along  the  curves  ^  = 
constant  and  0(  =  constant  respectively,  the  differential  forms 
E,  G  and  H  are  equal  to  unity.  This  simplifies  the  boundary 
conditions  and  considerably  reduces  the  numerical  effort  required 
to  obtain  a  solution  to  a  specific  problem. 

In  all  of  the  previous  development  It  was  assumed  that  the 
elastic  medltim  was  finite  and  bounded  by  a  closed  surface.  If 
the  medium  Is  Infinite  In  extent,  the  Integral  representations 
must  be  suitably  modified  and  the  displacement  potentials 
must  vanish  at  Infinity  In  a  prescribed  msmner. (9)> (12), (13) 

The  existence  of  solutions  for  three-dimensional  problems  Is 
discussed  by  Sokolnlkoff  and  Gurtln,  among  others,  who  give 
several  references . 


Discussion  and  Conclusions 


By  employing  methods  analogous  to  those  used  In  potential 
theory,  the  fundamental  boundary  value  problems  of  static  and 
dynamic  elasticity  have  been  formulated  In  terms  of  a  coupled 
set  of  Integral  equations  and  Inhomogeneous  partial  differential 
equations.  One  of  the  advantages  of  this  formulation  Is  that  the 
shape  of  the  boundary  Is  directly  accounted  for.  If  the  geometry 
of  the  cross-section  coincides  with  a  coordinate  system  In  which 
the  equations  satisfied  by  the  displacement  potentials  separate. 

It  Is  possible  to  obtain  solutions  by  eigenfunction  techniques. 
Another  advantage  Is  that  because  of  the  similar  structure  of 
the  formulation  of  each  problem  methods  developed  for  the  solution 
of  one  problem  may  more  easily  be  adapted  to  obtain  solutions  for 
the  others.  For  example.  If  a  numerical  method  of  solution  Is 
attempted,  the  Integrands  occurlng  In  all  of  the  formulations 
have  the  same  type  of  singularity.  Thus,  a  technique  that 
satisfactorily  estimates  the  contribution  to  the  Integral  In  the 
neighborhood  of  the  singularity  of  one  of  the  Integrals  may  be 
used  In  estimating  the  contributions  of  the  analogous  singularities 
In  the  other  Integrals. 

If  in  the  static  and  harmonic  time  dependent  problems  the 
Integrals  are  approximated  by  finite  sums,  the  resultant  set  of 
linear  equations  Is  simultaneous.  A  similar  numerical  approxima¬ 
tion  applied  to  the  Integrals  In  the  Initial  boundary  value  problem 
results  in  a  set  of  linear  equations  which  are  not  simultaneous. 
Hence,  the  implicit  and  explicit  cheu?acter  of  the  two  types  of 
problems  Is  clearly  evident.  Also,  It  should  be  noted  that,  the 
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first  amd  second  fundamental  boundary  value  problems  In  dynamic 
elasticity  are  analogous  to  the  problems  of  determining  the 
scattered  field  due  to  an  elastic  wave  Impinging  upon  an  arbitrary 
fixed  or  free  surface. 

After  the  desired  surface  potentials  are  determined,  the 
displacement  and  stress  fields  are  obtained  by  differentiating 
the  kernels  of  the  Integrals  and  properly  combining  the  resulting 
expressions  with  the  known  surface  potentials  and  Integrating. 

This  method  of  computing  the  stress  and  displacement  field  avoids 
the  decrease  In  rate  of  convergence  of  the  Infinite  series  that 
results  when  eigenfunction  expansions  are  used  to  obtain  these 
fields. 
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